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Abstract 

The  thermal izat ion  of  coherent  multiphoton  excitation  of  molecules  is 
examined.  An  inspection  of  the  level  distribution  of  a  driven  anharmonic 
oscillator,  which  is  subject  to  energy  and  phase  dissipation  to  a  bath,  reveals 
that  phase  relaxation  (dephasing)  due  to  the  bath  perturbation  play3  a  major 
role  in  the  thermal izat ion  process. 
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I.  Introduction 

Rapid  advances  in  laser  photochemistry  and  isotope  separation  have  generated 

considerable  interest  in  the  selective  absorption  of  electromagnetic  radiation 

i  2 

leading  to  either  molecular  dissociation  ’  or  adspecies  desorption  from  an 
3-5 

interface.  The  theoretical  treatment  of  the  problem  of  many  degrees  of  freedom 

for  a  polyatomic  molecule  and/or  admolecule-surface  system  is  very  complicated  and 

generally  an  insoluble  problem.  Besides  the  complexity  of  multiphoton  excitation 

which  is  nonresonant  for  high  transitions  due  to  the  presence  of  molecular 

anharmonicity,  interactions  between  the  resonant  photon-active  vibrational  mode  and 
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the  other  inactive  modes,  and/or  between  the  admolecule  and  surface  modes,  l.e., 

surface  phonons,  surface  plasmons,  reflected  photons  (photons  emitted  by  the 

admolecule  and  reflected  by  the  surface),  etc.,  create  a  large  number  of  dissipative 

channels  by  which  selective  absorption,  dissociation  and  desorption  at  the  surface 

are  strongly  influenced.  A  major  motivation  for  many  of  the  studies  in 

photochemistry  is  the  attractiveness  of  mode-selective  chemistry:  under  appropriate 

conditions,  molecules  prepared  in  a  state  of  high  vibrational  energy  by  absorption  of 

monochromatic  laser  radiation  would  be  markedly  different  from  the  ergodic  ensemble 

which  results  if  the  same  energy  is  imparted  by  heating,  i.e.,  by  collisions  with  a 

9.10 

thermal  bath.  Unfortunately,  experimental  and  theoretical  studies  have  shown  ’ 
that  near  the  dissociation  threshold  for  the  multiphoton  excitation  process  of  a 
polyatomic  molecule  or  near  the  desorption  threshold  for  such  excitation  of  an 

5 

admolecule-surface  system,  a  statistical  thermodynamic  process  occurs,  so  that  the 
dissociation  (or  adsorption)  reaction  proceeds  in  a  roughly  similar  manner  as  for  a 
thermally-excited  unimolecular  reaction.  In  this  letter,  we  use  a  simple  model  for  a 
laser-driven  damped  anharmonic  oscillator  for  the  active  vibrational  mode  which  is 
subject  to  both  energy  and  phase  dissipation  to  a  bath  (the  inter-  or  Intramolecular 
nonactive  modes  and  the  surface  excitations)  to  discuss  the  role  of  dissipative 
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mechanisms  in  the  thermal izat ion  coherent  multiphoton  excitation  processes.  We  shall 
see  that  phase  relaxation  (dephasing)  is  more  important  than  energy  relaxation  in  the 
active-mode  thermal izat ion  process.  The  transformation  of  the  level  population 
distribution  to  the  Planck  distribution  in  the  driven  damped  anharmonic  active  mode 
is  mainly  due  to  the  phase  dissipation  provided  by  the  bath. 

II.  The  Model 

The  excitation  and  relaxation  of  the  active  vibration  mode  in  a  collisionless 

multiphoton  process  can  be  visualized  in  terms  of  an  anharmonic  oscillator  which  13 

driven  by  a  near  resonant  coherent  laser  field  and  controlled  by  both  energy  (T^  and 

phase  (T2)  relaxation.  The  many  other  inactive  modes  and/or  surface  excitations 

(phonons,  plasmons,  reflected  photons,  etc.)  are  treated  as  a  large  thermal  bath,  by 

which  the  dissipative  mechanisms  are  provided.  A  reduced  density  operator  W  for  the 

12 

resonant  mode  obeys  the  master  equation 


dt  ”  ^"iLc+A1 +a2^w  • 


The  first  term  in  Eq.  (1),  -1L  W  ■  -( 1/fi) [H,W] ,  is  a  coherent  (reversible)  term, 

c  • 

where  H  -  H.  +  H_  is  the  sum  of  molecular  anharmonic  oscillator  and  driving  field 
M  r 


interactions. 


Hu  »  "fUw-w,  )a*a  -  fiea^a(a^a+1 ) 

M  L 


Hp  -  pE(t) (a+a) .  (3) 

Here  a*  and  a  are  oscillator  ladder  operators,  p  is  the  molecular  dipole,  E(t)  is  the 
applied  field,  u>L  is  the  laser  frequency,  and  e  is  the  (positive)  molecular 

anharmonicity.  The  last  two  terms  in  Eq.  (1)  are  dissipative  terms,  which  are  given 

.  11.12 

hir  * 


W  -  Y+([a,Wa*]  ♦  [aW,a*])  +  Yt([a*,Wa]  ♦  [a*W,a]) 


V:.:- 

.  '  V V  W  T  .  t  . 


y-.  .  •  »  , 


' -j" 


k 


A2W  -  ntCa^ajWa^al+Ca^aW, a*a])  .  (5) 

A2W  is  responsible  for  a  pure  phase  relaxation  mechanism  (dephasing),  Yt,  Y+  and  n 
contain  implicit  information  about  the  bath,  and  A^ W  is  responsible  for  the  energy 
relaxation  mechanism,  which  governs  the  population  loss  in  the  active  mode.  In  the 
Born  approximation,  i.e.,  where  the  characteristic  rate  of  the  coherent  exchange 
induced  by  the  system  Hamiltonian  is  sufficiently  smaller  than  the  largest  decay  rate 
of  the  dissipation  process,  the  diagonal  parts  of  the  density  matrix  obey 


d_ 

dt 


W 


n 


ft  -iL  (t~r)  A?(t-t) 

<n|A1Wp|n>  -  j  di  <nJPLIe  e  L^U)^  , 


(6) 


where  Lj  is  a  Liouville  operator  due  to  the  oscillator-field  interaction,  and  P  is  a 
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projection  operator  used  to  obtain  the  diagonal  part  of  the  density  operator. 

To  solve  the  difficult  problem  of  explicitly  evaluating  the  matrix  element  that 

appears  within  the  integrand  in  Eq.  6,  we  have  used  the  Weidlich  eigenfunction- 

expansion  method.12,114  The  "eigenstates”  of  the  super  operators  A2  and  A*  are  |  |Apq) 

and  | | Bpq ) ,  respectively.  For  a  normal  situation  in  molecular  excitation,  the  phase- 
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relaxation  rate  n  i3  always  much  larger  than  the  energy  relaxation,  i.e.,  n  >>  k, 
so  that  the  Markoff  approximation  is  applicable.  Within  this  approximation,  the 
equations  of  motion  for  the  diagonal  parts  of  the  density  matrix  are  then  given  as 


*  ft  _  2,  ♦  -iL  (  t~T ) 

-  <"U,«rl'»  -  )  dtlpq®  q  Tr(Bp(,LI"<*>><"lLie  *p,l">  • 


Equation  (7)  can  be  explicitly  written  as 


Vn  -  2«(i)[(n.1)Vln„-nVlnJ  •  2«(i  -1)[nWn.,-(n.1)Wn3 


<i  -u 


♦  r  (-1 )m+p(-i-)m  P^p"m*0)(22-T)/(p+1)]C(n,p)|  , 

‘•p-tD+l  1  "Z  D 


(8) 


where  G(n,p)  -  F(u-,  t)  (1  -z)"*1  nP^"1  1 1  }  (2z-1 )  -  F(uJ,t)  (1-z)n+2(n*1  )P^n  P,1)(2z-1)], 

F  (u^  »  t )  =  -  !"  2  n  /  (  2  )1  fe  nt(cos  w-t  -  (w*/n)sinu)-t)  -  ll 

n  L  nJI-  nn  nJ; 

w-  -  2( n-1 ) e  -  A,  w+  -  2ne  -  A,  z  -  (nB*1)  \  ng  is  the  averaged  excitation  of  the  bath, 
which  is  dependent  on  temperature,  A  is  the  laser  detuning  with  respect  to  the  fundamental 
frequency  of  the  oscillator,  P^  *  (x)  are  Jacobi  polynomials  of  the  indicative  argument, 

and  fiR  is  Rabi  frequency. 

For  the  case  of  low  bath  temperature,  i.e.,  where  the  average  excitation  of  bath 
quantum  nR<<1  and  the  rate  of  phase  relaxation  is  much  larger  than  the  energy  relaxation 
and  the  Rabi  frequency,  i.e.,  n  >>  <,  fiR,  Eq.  (8)  can  be  evaluated  by  a  lengthy  but 
straightforward  algebraic  calculation.  This  leads  to  the  following  set  of  coupled 
differential  equations,  which  takes  both  the  laser  coherent  excitation  and  the  thermal 
excitation  and  relaxation  into  account: 


Wn  '  2<(^+1)[(n*1)Wn+1-nWn]  ♦  Z^CnW^-Cn*!  )WR] 


?  n(W  -W  )  ( n+1 )  (W  -W  .) 

2nRp(t)  J— 5 - —  n  ^  )  . 


n2*[2e(n-1)-A]2 


n2+(2en-A)2 


(9) 


The  first  two  terms  of  Eq.  (9)  correspond  to  low- temperature  thermal izat ion 
processes  due  to  the  energy  relaxation,  which  causes  the  active  mode  to  approach  the 


steady  state  by  sharing  its  energy  with  other  degrees  of  freedom.  The  last  two  terms 
correspond  to  coherent  processes  due  to  the  laser  excitation,  where  the  influence  of 


mmc 
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the  bath  still  cannot  be  neglected  since  the  phase  relaxation  parameter  n  is  governed 
by  the  bath  interaction. 

1 5 

As  is  already  known,  the  level  population  distribution  of  a  pure  driven 

harmonic  oscillator  is  a  Poisson  distribution,  i.e.,  the  oscillator  has  a  single 

degree  of  freedom  with  no  coupling  to  the  other  degrees  of  freedom  and  its  initial 

state  is  coherent.  It  is  further  known  that  the  level  population  distribution  of  a 

single  driven  anharmonlc  oscillator  can  also  be  approximately  described  by  a  Poisson 

distribution.1^  As  for  a  fully  statistical  thermalized  system,  i.e.,  an  incoherently 

driven  quantum  oscillator,  the  level  population  should  be  in  a  Planck  distribution  at 

a  well-defined  temperature  determined  by  the  total  energy  deposited  in  the 
1 7 

oscillator.  It  is  hence  quite  interesting  to  know  how  the  population  distribution 
looks  for  an  anharmonic  vibrator  which  is  coherently  driven  by  a  laser  field  but  is 
subject  to  (incoherent)  energy  and  phase  dissipation. 

III.  Numerical  Results 

What  kind  of  dissipative  channel  is  more  important  in  the  thermalization 
processes  for  molecular  multiphoton  excitation?  To  answer  this  question,  we  evaluate 
the  vibrational  level  population  of  a  driven  damped  anharmonic  oscillator  subject  to 
both  energy  and  phase  relaxation  [Eqs.  (8)  and  (9)]  and  compare  it  with  the  Planck 
distribution  and  the  Poisson  distribution  having  th'  same  average  vibration 

—  OB 

excitation,  n  -  l  .  nW  .  We  see  that  when  the  phase  relaxation  rate  constant  n  is 
rather  small  (i.e.,  n*5  in  Fig.  1),  the  population  distribution  over  the  vibrational 
levels  of  the  driven  damped  anharmonic  oscillator  (D  curve)  is  markedly  different 
from  the  Poisson  distribution  (P  curve)  and  the  thermal  Planck  distribution  (T 
curve).  However,  if  we  let  the  phase  relaxation  constant  become  larger,  (e.g. ,  n-15 
in  Fig.  2),  the  curves  for  the  driven  damped  oscillator  (D  curves)  become  closer  to 
the  thermal  distribution  curves  (T)(see  Fig.  2).  We  also  notice  that  the  curves  for 
the  driven  damped  oscillator  never  behave  as  a  Poisson  distribution,  hinting  that  the 
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coherence  has  been  diminished  in  the  excitation  due  to  the  existence  of  the 
relaxation  pathway. 

To  determine  the  effect  of  energy  relaxation  on  the  population  distribution,  we 
have  also  computed  curves  for  different  values  of  the  energy  relaxaton  parameter 
t<»0. 1 ,G . 5,0 .9 , etc . ) ,  but  with  the  other  parameters  remaining  the  same  as  in  Fig.  1. 
Given  the  resulting  constant  proximity  of  the  D  curve  to  the  fully  thermalized  T 
curve,  we  conclude  that  energy  relaxation  has  a  minor  influence  on  the  thermalization 
of  the  driven  anharmonic  oscillator. 
l\l .  Summary 

We  have  used  a  model  for  a  driven  anharmonic  oscillator,  which  is  subject  to 
energy  and  phase  dissipation,  in  order  to  evaluate  thermalization  effects  in  the 
initial  stage  of  coherent  multiphoton  excitation.  It  is  shown  that  the  main  effect 
in  the  thermalization  of  the  vibrational  level  population  in  the  active  mode  is  due 
mainly  to  the  dephasing  agitation  coming  from  the  bath  interaction.  In  a  normal 
situation,  the  cephasing  rate  n  is  very  large  (n  1  is  on  the  order  of  a  picosecond), 
and  it  is  much  larger  than  the  energy  relaxation  rate  k.  To  achieve  mode-selective 
chemical  reactivity,  an  ultra-short  pulse,  i.e.,  sub-picosecond  laser  pulse,  may  be 
necessary.  The  mode-selective  reaction  must  take  place  in  the  period  of  the  short 
duration  of  laser  pulse  in  order  to  prevent  fast  dephasing. 
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Figure  Captions 

1.  Level  population  distribution  for:  (a)  a  laser-driven  damped  anharmonic 

oscillator  (D  curve);  (b)  Poisson  distribution  (P  curve);  (c)  thermal  Planck 

distribution  (T  curve).  Cases  (a)-(c)  all  have  a  same  averaged  excitation  energy 

n.  The  parameters  needed  to  calculate  the  D  curve  are  (k,  n,  0R,  e,  A)  -  (0.1, 

5,  1.0,  1.5,  3.0),  where  all  the  rates  are  in  the  unit  of  the  transverse 
1 2 

relaxation  rate.  n  must  be  a  positive  integer,  n  -  0,1,2,...,  and  there  are  no 
values  for  Wn  between  adjacent  n. 

2.  The  same  as  Fig.  1,  except  n- 1 5 . 
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